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Abstract. 

Interstellar dust grains exhibit complex dynamics which is essen- 
tial for understanding many key interstellar processes that involve dust, 
including grain alignment, grain growth, grain shattering etc. Grain ro- 
tational and translational motions are affected not only by gaseous col- 
lisions, but also by interactions with ions, photons, magnetic fields etc. 
Some of those interactions, e.g. interactions of ions with the dipole elec- 
tric moment of dust grains, require the quantum nature of the process to 
be accounted for. Similarly, coupling of rotational and vibrational degrees 
of freedom in a grain happens due to relaxation processes, among which 
the process related to nuclear spins frequently is the dominant one. This 
coupling modifies substantially both the dynamics of rotational and trans- 
lational motions by inducing grain flips. The flipping averages our certain 
systematic torques and forces that act on grains. As the rate of flipping 
is larger for smaller grains, these grains can be "thermally trapped", i.e. 
rotate at a thermal rate in spite of the presence of systematic torques. 
Moreover, a subset of small grains with high dipole moments may rotate 
at subthermal rates due to high damping arising from grain emission in 
microwave range of frequencies. Translational and rotational dynamics of 
grains is interrelated. For instance, both rotation and gas grain relative 
speeds affect grain alignment as well as determine grain size distribution 
and structure. Translational dynamics of grains is mostly dominated by 
grain interactions with magnetohydro dynamic turbulence. Efficient tur- 
bulent mixing of dust grains limits the degree to which grains of different 
sizes may be segregated in space. 



1. Why do we care? 

Dust is an important constituent of the interstellar medium (ISM). It interferes 
with observations in the optical range, but provides an insight to star-formation 
activity through far-infrared radiation. As dust gets aligned in the external 
magnetic field (Hall 1949; Hiltner 1949) it traces the magnetic field via emission 
and extinction polarization (see reviews by Hildebrand et al. 2000; Lazarian 
2003). Grain alignment is a property of rotating grains. Fast rotation makes 
grains less susceptible to disorientation by collisions with gas atoms, but it also 
limits the efficiency of some alignment mechanisms. In addition, fast rotating 
grains are the sources of microwave emission that interferes with measurements 
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of Cosmic Microwave Background (CMB) intensity and polarization ( Draine &i 
Lazarian 1998a; see also recent review by Lazarian &; Finkbeiner 2003). 

The basic properties of dust (optical extinction, dust chemistry, heating of 
the ISM etc) strongly depend on its size distribution (Biermann & Harwit 1980; 
O'donnell &: Mathis 1997). The latter evolves As a result of grain collisions, 
whose frequency and consequences depend on grain relative velocities (Draine 



Rotational and translational motions are not independent. For instance, the 
difference in the number of active sites of H2 formation over the parts of the grain 
perpendicular to the axis of grain rotation can result in uncompensated force 
accelerating the grain (Purcell 1979, Lazarian & Yan 2002, henceforth LY02). 
However, if a grain is undergoing frequent flips due to thermal fluctuations 
(Lazarian Sz Draine 1999a, henceforth LD99a and §2) the thrust will change 
direction and be averaged out. On the other hand, grain alignment presents a 
case when translational motion affects grain rotation. For instance, mechanical 
alignment of grains (see Lazarian 2003 and references therein) induces higher 
rotational rates if grains rotate thermally initially. 

Damping of rotational and translational motions have many similarities. 
For instance, the interaction of grains having dipole moments with passing ions 
is important for damping of both types of motions. For the case of damping of 
translational motion the treatment of the process is given in Yan, Lazarian & 
Draine (2003, henceforth YLD03). Quantum cut-off discussed in §4 is applicable 
to this case as well. 

In what follows we discuss how dust grains rotate (§2), at what rate they ro- 
tate (§3), and whether quantum effects are important for grain rotation (§4). We 
describe the translational motion of grains in §5. In §6 we discuss astrophysical 
implications of grain dynamics and provide the summary in §7. 

2. Do grains rotate about their axis of major inertia? 

Originally this question was asked in relation with the theory of grain alignment 
(see Fig. la and discussion in Lazarian 2003). To produce the observed starlight 
polarization, grains must be aligned with their long axes perpendicular to mag- 
netic field. This involves alignment not only of grain angular momentum J in 
respect to the external magnetic field B, but also the alignment of grain long 
axes in respect to J. Jones Sz Spitzer (1967) assumed a Maxwellian distribution 
of angular momentum which favored the preferential alignment of J with the 
axis of the maximal moment of inertia (henceforth axis of major inertia). Pur- 
cell (1979, henceforth P79) later considered grains rotating much faster than the 
thermal velocities (sec §3) and showed internal dissipation of energy in a grain 
will make grains rotate about the axis of major inertia. 

Arguments in P79 can be easily understood. Indeed, for an oblate grain 
(see Fig. lb) with angular momentum J the energy can be written as 



where h = Imax/Ii. is the ratio of the maximal to minimal moments of grain 
inertia. Internal forces cannot change the angular momentum, but it is evident 
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from eq.(l) that the energy can be decreased by ahgning the axis of maximal 
inertia along J, i.e. by decreasing (3. P79 discussed two possible causes of 
internal dissipation, the first one related to the well known inelastic relaxation 
(see also Lazarian &: Efroimsky 1999), the second is due to the mechanism that 
he discovered and termed "Barnett relaxation". 

We remind the reader that the Barnett effect is converse of the Einstein-de 
Haas effect. If in Einstein-de Haas effect a paramagnetic body starts rotating 
during remagnetizations as its flipping electrons transfer the angular momentum 
(associated with their spins) to the lattice, in the Barnett effect the rotating body 
shares its angular momentum with the electron subsystem causing magnetiza- 
tion. The magnetization is directed along the grain angular velocity and the 
value of the Barnett-induced magnetic moment is /i ~ 10^^^a;(5) erg gauss"^ 
(where a;(5) = a;(s"^)/10^)^ 

Simplified Model of Alignment 



Figure 1. Left panel- Grain alignment implies several alignment pro- 
cesses acting simultaneously and spanning many time scales (shown for 
10~^ cm grain in cold interstellar gas). The rotational dynamics of a 
grain is rather complex. The internal alignment introduced by Pur- 
cell (1979) was thought to be slower than precession until Lazarian & 
Draine (1999b, henceforth LD99b) showed that it happens 10^ times 
faster when relaxation through induced by nuclear spins is accounted 
for (approximately 10^ s for the 10~^ cm grains). Right panel- Grain ro- 
tation arising from systematic torques arising from H2 formation (P79). 
In the presence of efficient internal relaxation the angle f3 between the 
axis of maximal moment of inertia and J is small is grain is rotating 
at suprathermal rates (E'fcmetic > kTgrain)- 



^Therefore the Larmor precession has a period ~ 3 x lO^B^gj s (where ^) and the magnetic 
field defines the axis of ahgnment (see more details in Lazarian 2003) 




Precession of J is rapid => magnetic field 
is the axis of alignment 
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Figure 2. Left panel- Schematic of a grain affected by Purcell's 
torques before and after a flipover event. As the grain flips the di- 
rection of torques changes. H2 formation sites act as thrustcrs in the 
figure. Right panel- Regular crossover event as it described by Spitzer 
&; McGlynn (1979). The systematic torques decrease the amplitude 
of J component parallel to the axis of maximal inertia to zero while 
preserving the other component, J±. If J± is small then the grain is 
susceptible to randomization during crossovers. The direction of J is 
preserved in the absence of random bombardment. 



The Barnett relaxation process may be easily understood. We know that a 
freely rotating grain preserves the direction of J, while angular velocity precesses 
about J and in grain body axes. The "Barnett equivalent magnetic field", i.e. 
the equivalent external magnetic field that would cause the same magnetization 
of the grain material, is Hbe = 5.6 x 10~^a;(5) G. Due to the precession of 
angular velocity the "Barnett equivalent magnetic field" changes in grain axes. 
This causes remagnetization accompanied by the inevitable dissipation. As a 
result E{9) and therefore decreases. 

The Barnett relaxation happens on the scale tsar ~ 4 x IO'^lu'^^ sec, very 
short compared to the time tgas over which randomization through gas-grain 
collisions takes place. As a result, models of interstellar polarization developed 
since 1979 have often assumed that Barnett dissipation aligns J perfectly with 
the major axis of inertia. However, Lazarian (1994, henceforth L94) has pointed 
out that this is a poor approximation if the grains rotate with thermal kinetic 
energies: thermal fluctuations in the Barnett magnetization will excite rotation 
about all 3 of the body axes, preventing perfect alignment unless the rotation 
is suprathermal or the grain solid temperature is zero. The exact analysis of 
the problem was given in Lazarian &i Roberge (1997), where the distribution 
of 6 for a freely rotating grain was defined through the Boltzmann distribution 
exp(—E(9)/kTgrain), where Tgrain is the grain temperature. 
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3. How fast do the grains rotate? 

Earlier work in the field assumed Brownian grain rotation with the eff^ective 
temperature equal to the mean of the grain and gas temperatures (see Jones & 
Spitzer 1967). The complexity of grain rotation was realized only later. Purcell 
(1975; 1979) realized that grains may rotate at a much faster rate resulting from 
systematic torques. P79 identified three separate systematic torque mechanisms: 
inelastic scattering of impinging atoms when gas and grain temperatures differ, 
photoelectric emission, and H2 formation on grain surfaces (see Fig. lb); we will 
refer to these below as "Purcell's torques". The latter was shown to dominate 
the other two for typical conditions in the diffuse ISM (P79). The existence of 
systematic H2 torques is expected due to the random distribution over the grain 
surface of catalytic sites of H2 formation, since each active site acts as a minute 
thruster emitting newly-formed H2 molecules. 

Independent of Purcell, Dolginov (1972) and Dolginov & Mytrophanov 
(1975) identified radiative torques as the way of spinning up a subset of interstel- 
lar grains. "Helical" grains would scatter differently left and right polarized light 
and therefore ordinary unpolarized light would spun them up. This subset of 
"helical" grains was believed to be somewhat limited to special shapes/materials. 
This work did not make much impact to the field until Draine & Weingartner 
(1996) showed that grains of arbitrary irregular shapes get spun up when their 
size is of the order of the starlight wavelength. 

The arguments of P79 in favor of suprathermal rotation were so clear and 
compelling that other researchers were immediately convinced that interstellar 
grains in diffuse clouds should rotate suprathermally. Purcell's discovery was of 
immense importance for grain alignment. Suprathermally rotating grains remain 
subject to gradual alignment by paramagnetic dissipation (Davis &; Greenstein 
1951), but due to their large angular momentum are essentially immune to 
disalignment by collisions with gas atoms. 

Spitzer Sz McGlynn (1979, henceforth SM79) showed that suprathermally 
rotating grains should be susceptible to disalignment only during short inter- 
vals of slow rotation that they called "crossovers" (see Fig 2, right). During 
a crossover the grain slows down, fiips and is accelerated again. Crossovers 
are due to various grain surface processes that change the direction (in body- 
coordinates) of the systematic torques. Therefore for sufficiently infrequent 
crossovers suprathermally rotating grains will be well aligned with the degree of 
alignment determined by the time between crossovers, the degree of correlation 
of the direction of grain angular momentum before and after a crossover (SM79) , 
and environmental conditions (e.g., magnetic field strength B). 

The original calculations of SM79 obtained only marginal correlation of 
angular momentum before and after a crossover, but their analysis disregarded 
thermal fluctuations within the grain material. Indeed, if the alignment of J with 
the axis of major inertia is perfect all the time through the crossover the absolute 
value of |J| passes through zero during the crossover. Therefore gas collisions 
and recoils from nascent H2 molecules would completely randomize the initial 
and final directions of J during the crossover. Thermal fluctuations partially 
decouple J and the axis of major inertia (see §2). As a result the minimal value 
of |J| during a crossover is equal to the component of J perpendicular to the 
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Figure 3. Left panel- Grain trajectory on the J± - J\\ plane, where 
Jj_ and J|| are components of J perpendicular or parallel to the grain's 
principal axis of largest moment of inertia. The solid trajectory shows 
a "thermal flip", while the broken line shows the "regular" crossover 
which would occur in the absence of a thermal flip. From LD99a. 
Right panel- Top: Thermal flipping to damping ratio as a function 
of J/Jthermai for grains of given size [broken lines, labeled by a_5 = 
a(cm)/10~^] and for grains with a given value of systematic torques 
[solid lines, labeled by G]. Dot shows Jmin = J'^tf for flipping-assisted 
crossover of a_5 = 0.5 grain with G = 10^. Bottom: Thermal trapping 
for grains of given size [broken lines, labeled by a-s], and given value 
of torques [solid lines, labeled by G]. From LD99a. 

axis of major inertia. This value is approximately Jthermai ~ {'^kTgrainlm.axY^'^ 
and the randomization during a crossover decreases (Lazarian & Draine 1997, 
henceforth LD97). LD97 obtained a high degree of correlation of angular mo- 
mentum direction before and after the crossover for grains larger than a critical 
radius ~ 1.5 x 10~^cm, the radius for which the time for internal dissipation 
of rotational kinetic energy is equal to the duration of a crossover. 

What would happen for grains that are smaller than Oc? Lazarian &; Draine 
(1999a, henceforth LD99a) showed that instead of following the phase space 
trajectory prescribed by SM79 theory along which J_|_ is approximately constant, 
while the component of J parallel to the axis of maximal inertia Jy changes sign, 
the grains undergo flipovers (see Fig2a) during which | J| does not change (see 
FigSa). If these flipovers repeat; grains get "thermally trapped" (LD99a and 
FigSb). This process can be understood in the following way. For sufficiently 
small |J| the rate of flipping t^^ becomes large. Purcell's torques change sign 
as grain flips. As a result, grain undergoes a random walk over the time of 
grain rotational damping td- If Jth is the momentum of a thermally rotating 
grain, the expected value of the mean squared angular momentum (J^) will be 
approximately J^f^ + {G — l)J^f^ttf/{td+ttf), where G is defined by the systematic 
Purcell's torque being (G — 1)^^"^ Jth/td- As ttf{J'^) is a non-linear function of 
angular momentum the solution for the angular momentum is bistable. For 
sufficiently small initial less than the critical value J^^, it tends to « Jj^, 
while for of > Jf^ it tends to its suprathermal value predicted by Purcell 
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i.e. ~ GJfy^. A siiprathermally rotating grain has chances to become thermally 
trapped and stay in the trapped state for ~ td^^vii Jcr/ Jth)'^\ till fluctuations in 
the value of the angular momentum rescue the grain from the trap. As a result a 
substantial part of grains smaller than Ucr will not rotate at high rates predicted 
by P79 even in spite of the presence of systematic torques that are fixed in body 
coordinates. Thermally rotating grains are subject to randomization by gaseous 
collisions and arc marginally aligned for typical interstellar fields (see Lazarian 
1997a; Robcrge &: Lazarian 1999). Other consequences of thermal trapping and 
the values of Ucr that follow from processes different from the Barnett relaxation 
are discussed below. The LD99a model gives plausible predictions, but a more 
comprehensive quantitative study is necessary. We are aware of such a study 
being done by Roberge &: Ford (in preparation; see also the contribution by 
Roberge (2004) to this volume). 

Thermal trapping limits the range of grain sizes which can be spun up 
by Purcell's torques. One may expect that radiative torques are not much 
influenced by thermal trapping as they are not fixed in the grain axes and may 
not alter their direction due to grain flipping. A quantitative study by Draine & 
Weingartner (1996) shows that for typical interstellar spectra grains with sizes 
larger than 5 x 10~^ cm can be spun up by radiative torques^ This means that in 
ISM grains larger than 5 X 10~^ cm rotate suprathermally^, while small grains 
may rotate thermally. 

Can grains rotate subthermally, i.e. with velocities which are substantially 
smaller than those of Brownian rotation? The answer to this is positive. Draine 
k. Lazarian (1998a, henceforth DL98) showed that grains with sizes less than 
10~^ cm may be efficient emitters of microwave radiation provided that they have 
dipole moments. The emission from those grains can account for the so-called 
Foreground X (De Oliveira-Costa et al. 2002) observed in the range 10-100 GHz 
(see Fig. 4a). The back-reaction from radiation slows down grain rotation. As 
a result, detailed calculations in DL98 showed that emitting grains may rotate 
at substantially subthermal velocities, resulting in considerably less microwave 
emission than had been estimated by Ferrara & Dettmar (1994), who assumed 
Brownian rotation. 

4. Are quantum effects important for grain dynamics? 

Dust grains are essentially macroscopic bodies. The rotational quantum number 
J !%> 100 even for the smallest grains (DL98). Therefore the initial reaction is 
that grain dynamics can always be described using classical physics. This is not 
true, however. 

DL98 showed that to describe the interaction of a grain with a dipole mo- 
ment with ions one must account for the quantum nature of the interaction. 
Indeed, the pioneering work by Anderson & Watson (1993) assumed that ions 
within the Debye screening length interact with the grain dipole moment. 



^In the vicinity of stars with UV excess smaller grains can be spun up as well. 

^An important point here is that radiative torques may uot be dominant in terms of their 
absolute value, but still they can rescue grains from thermal trapping. 
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Calculations in DL98 showed that the more severe limitations on which ions 
can interact with grain dipole moment come in case of sufficiently small grains 
from the requirement that the angular momentum transfered in an individual^ 
interaction should be multiples of h. 

Another case when quantum mechanics is important is related to the Bar- 
nett effect. The Barnett effect is a quantum effect, as it involves orientation 
of electron spins, which are quantum objects. P79 noted that an analog of the 
Barnett effect exists for nuclear spins. If a rotating body has initially an equal 
number of nuclear spins directed parallel and anti-parallel to the angular veloc- 
ity U, it can decrease its kinetic energy, at constant total angular momentum 
J, if some of the angular momentum is transfered to the nuclear spin system. 
Increasing the projection of the nuclear angular momentum along J by +h (at 
constant J) reduces the rotational kinetic energy by HQ. If the rotating body is 
allowed to come into thermal equilibrium (without exchanging angular momen- 
tum) with a heat reservoir of temperature T^ust then particles of spin S develop 
a net alignment per particle 

Em=-5 e^pjmhn/kTdust) 



■im= 



_gex.p{m.hn/kTdust) 



Note that this docs not depend on the magnetic moment /i. 

As the number of parallel and antiparallel spins becomes different the 
body develops magnetization. The relation between and the strength of the 

(n) 

"Barnett-equivalent" magnetic field H^^^ (P79) that would cause the same nu- 
clear magnetization (in a nonrotating body) is given by 

where gn is the so-called nuclear ^-factor (see Morrish 1980), and /xn = eh/2mpC 
is the nuclear magneton, smaller than the Bohr magneton by the electron to 
proton mass ratio, nie/mp. A striking feature of eq. (3) is that the Barnett- 
equivalent magnetic field is inversely proportional to the species magnetic mo- 
ment. As grain tumbles the magnetization changes in grain body coordinates 
and this causes paramagnetic relaxation. This relaxation is proportional to 
X^iojJH]^^ (where x'n imaginary part of the nuclear contribution to the 

susceptibility) and is approximately 10^ times faster for nuclear moments than 
for their electron counterparts^. 

The first striking question is why an effect so feeble as nuclear magnetism 
can be so important in terms of internal relaxation. The answer is that spins 



^Excitation of plasma waves as an additional source of damping was considered in Ragot (2002) . 
However, it is possible to show that for small interstellar grains considered in DL98 collective 
plasma effects are negligible. 

® Obviously enough this sort of arguments fails when the rate of rotation is larger than the rate 
at which rmclcar spins interact. The latter depends on the concentration of nuclear spins. For 
the ISM dust LD99b concluded that the nuclear relaxation is dominant for thermally rotating 
grains larger than ~ 5 x 10~^ cm. 
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Figure 4. Left "panel- Emission of spinning grains for different phases 
of interstellar media in the model by Draine & Lazarian (1998b) against 
observational data (WMAP data is by D. Finkbeiner). From Lazarian 
& Finkbeiner 2003. Right panel - Time r for alignment of a with J for 
a freely rotating grain, due to dissipation associated with nuclear spins 
(nuclear relaxation) and electron spins (Barnett effect). Results are 
shown for grains with (J/J^)^ = 1 (bottom) and 10 (top). Also shown 
is the "crossover time" tc = J / where the torque J\\ is due to H2 
formation with a density of active sites 10^^cm~'^cm~^ (From LD99b). 



align in a rotating body because of their angular momentum, not because of 
their magnetic moments. The magnetic moment enters only as a means for the 
spins to exchange angular momentum with the lattice. The coupling within the 
electron spin system is very effective, with the result that there is minimal "lag" 
in the electron spin alignment when the grain angular velocity $7 changes in grain 
body coordinates. In the case of the nuclear spin system, however, the value of 
the nuclear magnetic moment is large enough to provide significant coupling, but 
weak enough so that there is a significant lag^ in the nuclear spin alignment when 
r2 precesses around the grain axis of maximal inertia: the coupling is "just right" 
for nuclear relaxation to be extremely effective for 10~^cm ^ a ^ 10~^cm grains. 
The efficiency of nuclear relaxation drops for sufficiently high frequencies and 
therefore Barnett relaxation dominates for a ^ 5 x 10~^cm, as seen in Fig. (4a). 

As we discussed in §3, LD99a found that "thermal flipping" was an impor- 
tant element of grain dynamics for sufficiently small grains. LD99b have shown 
that internal relaxation associated with nuclear spin alignment can be many 
orders of magnitude more rapid than due to the Barnett effect for a ^ 10~^cm 
grains. As a result, the phenomena of thermal flipping and thermal trapping 
become important for grains as large as ~ 10~^cm. 



''Purcell (1969) was the first to notice that for thermally rotating 10~^ cm grains the rate of 
paramagnetic alignment in the external magnetic field does not depend on the value of magnetic 
moment of the paramagnetic species. Therefore he concluded that the alignment of grains by 
interstellar magnetic field may happen due to nuclear spins only. 
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The Barnett effect is related to yet another quantum process relevant to 
the paramagnetic alignment of ultra-small grains (< 10~^ cm) that is discussed 
in Lazarian & Draine (2000). Those grains produce important CMB foreground 
(sec Fig. 4) which can be polarized if the grains are aligned. 

Introduced by Davis & Greenstein (1951) paramagnetic relaxation is easy 
to understand: for a spinning grain the component of interstellar magnetic field 
perpendicular to the grain angular velocity varies in grain coordinate. The 
resulting time-dependent magnetization has associated energy dissipation and 
torques that damp grain rotation perpendicular to magnetic field. 

Rotation removes the spin degeneracy of the electron energy levels. The 
energy difference between electron spin parallel or antiparallel to Q provides a 
level splitting corresponding to Ulo = ff/XB-f^BE- 

Now consider a (weak) static magnetic field H at an angle 9 to Cl. In grain 
coordinates, this appears like a static field HcosO plus a field HsinO rotating 
with frequency lo. This rotating field can be resonantly absorbed, since the 
energy level splitting is exactly huj. 

This energy absorption, of course, is proportional to paramagnetic suscepti- 
bility x"(a;). In the classical Davis-Greenstein analysis this magnetic susceptibil- 
ity x(u;) is taken to be that of a sample at rest in zero magnetic field: x ^ xi^o = 
0,u). Here we point out that one should instead use % = x±{Ho = Hbe,oj), 
where x±iHo,uj) describes the response of nonrotating material to a weak field 
rotating at frequency lo perpendicular to a static magnetic field Hq. Compared 
to classical Davis-Greenstein relaxation the study by Lazarian &; Draine (2000) 
predicts that the relaxation happens at the maximal possible or resonance rate. 
Therefore the relaxation was termed "resonance relaxation" . For grains of the 
order of 10~^ cm this results in a correction factor of order unity. For ultra-small 
grains it provides an enhancement of many orders of magnitude. 

5. Is turbulence important for driving grains? 

Depending on grain relative velocities, grain-grain collisions can have various 
outcomes, e.g., coagulation, cratering, shattering, vaporization and erosion and 
ejection of mantle material (see Draine 1985 and references therein). It is likely 
that some features of grain size distribution (Mathis, Rumpl & Nordsicck 1977; 
Kim, Martin & Hendry 1994), eg., cutoff at large size, are the result of fragmen- 
tation. 

Dust dynamics is also influential to the mctallicity in various astrophysical 
environment. The depletion of heavy elements in interstellar medium is a long 
puzzling problem. Grains moving supersonically may efficiently accrete gas- 
phase heavy elements (Weingartner Sz Draine 1999, Wakker Sz Mathis 2000). 
Dust ejection from galaxies is important to metal enrichment of intergalactic 
medium (see Aguirre et al. 2001 and references therein). 

Various processes can affect the velocities of dust grains. Shocks, radia- 
tion, ambipolar diffusion, and gravitational sedimentation all can bring about a 
dispersion in grain velocities (Draine 1985). 

Shock acceleration. The basic idea is that the weakly charged grains are 
like ions with high mass to charge ratio (Epstein 1980). Thus they can easily 
diffuse farther back upstream of the shock and be accelerated more efficiently to 
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suprathermal energies. Nevertheless, the shock acceleration is inefficient for low 
speed grains. The reason is that the efficiency of the shock acceleration depends 
on the scattering rate, which is determined by the stochastic interaction with 
turbulence. For low speed particles, the scattering rate is lower than the rate 
of momentum diffusion. In this case the stochastic acceleration by turbulence 
happens faster than dust acceleration by shocks (Park & Petrosian 1996, Yan Sz 
Lazarian 2003, henceforth YL03, see also below). 

Effects of radiation fields. Grains are exposed to various forces in 

anisotropic radiation fields. Apart from radiation pressure, grains are subjected 
to forces due to the asymmetric photon-stimulated ejection of particles. A de- 
tailed discussion can be found in Weingartncr & Draine (2001a, henceforth 
WDOl; also the review by Draine 2003). The photoelectric force depends on 
the ambient conditions relevant to grain charging (WDOl; Lafon 1990; Kerker 
&; Wang 1982). The calculation by WDOl demonstrated that it is comparable 
to radiation pressure when the grain potential is low and the radiation spec- 
trum is hard. Photodesorption is also a photon-stimulated ejection process, but 
of absorbed atoms on grain surface. However, because the surface physics and 
chemistry of grain materials are unclear, the calculation of the photondesorp- 
tion force is highly uncertain (WDOl). The bottom line is that the force due 
to photodesorption is expected to be comparable to the radiation pressure and 
photoelectric thrust (Draine 2003). 

Motions arising from H2 thrust. A different residual imbalance arises 
from the difference of the number of catalytic active sites for H2 formation on 
opposite grain surfaces. The nascent H2 molecules leave the active sites with ki- 
netic energy E, and the grain experiences a push in the opposite direction. The 
uncompensated force above is parallel to the rotational the other compo- 

nents of force are averaged out due to grain fast rotation. Applying the charac- 
teristic values^ in Lazarian & Draine (1997), LY02 got the "optimistic" velocity 
V ~ 330(10"^cm/a)i/2cin/s for CNM (Fig. 5a) and v ~ 370(10-5cm/a)°-^cm/s 
for WNM. The grains tend to be aligned with rotational axes parallel to the 
magnetic field. Thus grains acquire velocities along the magnetic field lines. 
It is clear from Fig. 5a that for the chosen set of parameters the effect of H2 
thrust is limited. The percentage of atomic hydrogen is reduced in dark clouds, 
and the radiation field is weak. Therefore, the velocities driven by the variation 
of the accommodation coefficient are always much smaller than those due to 
turbulence so that grains in dark clouds should be fully mixed. 

The interstellar medium is turbulent (see Arons & Max 1975; Scalo 1987; 
Lazarian 1999). Turbulence has been invoked by a number of authors (see 
Kusaka et al. 1970; Volk et al. 1980; Draine 1985; Ossenkopf 1993; Weiden- 
schilling &; Ruzmaikina 1994) to provide substantial grain relative motions. 

Grain Motions due to Frictional Drag 
In hydrodynamic turbulence, the grain motions are caused by the frictional in- 
teraction with the gas. At large scale grains are coupled with the ambient gas, 
and the slowing fluctuating gas motions mostly cause an overall advection of 
the grains with the gas (Draine 1985). At small scales grains are decoupled. 



^Tho number of H2 formation sites is highly uncertain. It may also depend on the interplay of 
the processes of photodesorption and poisoning (Lazarian 1995b; 1996). 



12 



Lazarian &; Yan 



The largest velocity difference occurs on the largest scale where grains are still 
decoupled. Thus the characteristic velocity of a grain in respect to the gas 
corresponds to the velocity dispersion of the turbulence on the time scale tdrag 
(Draine & Salpeter 1979). Using the Kolmogorov scaling relation oc k"^/^, 
Draine (1985) obtained the largest velocity dispersion in hydrodynamic turbu- 
lence V ~ V (tdrag /Tmax)^^^ , where Tmax is the time scale of the turbulence at 
the injection scale. 

As most interstellar arc magnetized and magnctohydrodynamic (MHD) tur- 
bulence is much better understood now, a revisit on the problem, namely accel- 
eration by MHD turbulence, is needed. In the following, we shall first introduce 
some basic idea of MHD turbulence and then give a brief review of our work 
on dust dynamics in MHD turbulence. MHD perturbations can be decomposed 
into Alfvenic, slow and fast modes (see Alfven & Falthmmar 1963). Alfvenic 
turbulence is considered by many authors as the default model of interstellar 
turbulence. Unlike hydrodynamic turbulence, Alfvenic turbulence is anisotropic, 
with eddies elongated along the magnetic field. This happens because it is easier 
to mix the magnetic field lines perpendicular to the direction of the magnetic 
field rather than to bend them. As eddies mix the magnetic field lines at the 
rate k±Vk , where Vk is the mixing velocity at this scale, the magnetic pertur- 
bations (waves) propagate along the magnetic field lines at the rate A;||Va • The 
Alfven and slow modes can be described by GS95 model (Goldreich & Sridhar 
1995; Cho &: Lazarian 2002; Lithwick &; Goldreich 2001; see also Cho, Lazarian 
Sz Vishniac 2002 for a review) . The corner stone of the GS95 model is a critical 
balance between these rates, i.e., k_iVk ~ k^\VA, which may be also viewed as 
coupling of eddies perpendicular to the magnetic field and wave-like motions 
parallel to the magnetic field. Calculations by Cho Sz Lazarian (2002; 2003c) 
demonstrated that fast modes arc very similar to acoustic turbulence. 

In the MHD case, grain motions are affected by magnetic fields. The charged 
grains are subjected to the electromagnetic forces which depend on the grain 
charge. If the periods of Larmor motion tl are longer than the gas drag time 
idragj the grains do not feel magnetic field directly. Otherwise, grain perpendic- 
ular motions are constrained by magnetic field. 

As Alfvenic turbulence is anisotropic, it is convenient to consider separately 
grain motions parallel and perpendicular to the magnetic field. The correspond- 
ing discussions can be found in LY02, Lazarian Sz Yan 2002b. 

The velocity dispersion induced by the compressional motion associated 
with the fast modes also causes the relative movement of the grain to the ambient 
gas (see YL03). The velocity fluctuations associated with fast modes are always 
in the direction perpendicular to B in a low (3 ((3 = Pgas/Pmag) medium. Thus 
the grain velocities are also perpendicular to B. In high (3 medium, grains can 
have velocity dispersion in any direction as the velocity dispersions of fast modes 
are radial, i.e., along k. 

Acceleration of Grains by Gyroresonance 
In YL03 we identified a new mechanism of grain acceleration, namely, gyrores- 
onance that is based on the direct interaction of charged grains with MHD 
turbulence. There exists an important analogy between dynamics of charged 
grains and dynamics of cosmic rays (see Yan Sz Lazarian 2002), and the ex- 
isting machinery developed for cosmic rays was modified to describe charged 
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Figure 5. Left panel- Relative velocities as a function of radii in 
CNM for silicate grains. The dotted line represents the gyroresonance 
with fast modes. The dashdot line refers to the gyroresonance with 
Alfven modes. The cutoff is by viscous damping. The dashed line is 
the result from hydro drag with Alfven modes (sec LY02), the solid line 
represents the hydro drag with fast modes. Contributions from different 
mechanisms are approximately additive in squares. The grain velocity 
driven by H2 formation (solid line). The part marked by open circles is 
nonphysical because thermal flipping is not taken into account (YL03). 
Right panel- Grain velocities in CNM gained from gyroresonance for 
different magnetic field strengths (YLD03). 

grain dynamics. The energy exchange involves resonant interactions between 
the grains and the waves (see YL03). Specifically, the resonance condition is 
LO — k^^vji = nil, (n = 0, ±1, ±2...), which means that the Doppler-shifted fre- 
quency of the wave in the grain's guiding center rest frame Ugc = 00 — k\\viJ, is 
a multiple of the particle gyrofrcqucncy Q. Basically there arc two main types 
of resonant interactions: gyroresonance acceleration and transit acceleration. 
Transit acceleration (n = 0) requires longitudinal motions and only operates 
with compressible modes. As the dispersion relation is w = kVf > kV^ for fast 
waves, it is clear that it can be only applicable to super-Alfvenic (for a low /? 
medium) or supersonic (for a high /? medium) grains. For low speed grains that 
we deal with here, gyroresonance is the dominant MHD interaction. 

The calculation by YL03 showed grains gain the maximum velocities per- 
pendicular to the magnetic field and therefore the averaged decreases. This 
is understandable since the electric field which accelerates the grain is in the 
direction perpendicular to the magnetic field. 

Then we applied our results to various idealized phases of interstellar medium 
(YLD03) . In Fig. 5a, we show the velocity of grain as a function of grain size in 
CNM. 

The acceleration by gyroresonance in both MC and DC arc not so efficient as 
in other media. This happens because the time for the gyrorcsonant acceleration 
tdrag a-^e much shorter in MC and DC. In these media, the drag time because of 
the high density is less than the gyro-period for grains larger than 10~^cm. 

It should be noted that the strength of magnetic fields in ISM is still some- 
what uncertain and may vary from place to place. We adopted a particular set 
of values in above calculations. How would the results vary as the magnetic 
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field strength vary? First of all, we know that the critical condition for acceler- 
ation is that the Larmor period tl is greater than the period of the turbulent 
motions at the damping scale. Grains with tl less than that can not be acceler- 
ated via gyroresonance. Thus the cutoff grain mass rric varies with the medium, 
iTT-c ~ QTcB I c. The magnitude of the velocity is a complex function of the mag- 
netic field. For illustration, here we demonstrate the results for 3 times stronger 
and weaker magnetic field (see Fig. 5b). Since the acceleration via hydro drag 
by fast modes decreases with the magnetic field, the relative importance of the 
two processes, namely, gyroresonance and hydro drag depends on the magnitude 
of the magnetic field. In magnetically dominant region, gyroresonance is domi- 
nant. In weakly magnetized region, the decoupled motions provide the highest 
acceleration rate. The injection scale is another uncertain parameter, but the 
grain velocity is not so sensitive to it provided that the injection scale is much 
larger than the damping scale. 

6. Discussion 

Implications of Rotational Dynamics 

Grain Alignment 

Grain alignment critically depends on rotational dynamics of grains. In- 
corporation of thermal fluctuations within theories of grain alignment (Lazarian 
1994, 1997ab, Lazarian & Roberge 1997, Lazarian & Drainc 1997, Roberge &: 
Lazarian 1999, Weingartner &; Draine 2003) introduced substantial changes in 
our understanding how grains get aligned. Thermal flips (LD99a,b) and thermal 
trapping (LD99a) that they entail have radically changed some of the contem- 
porary paradigms of grain alignment. Indeed, if grains get thermally trapped 
as discussed earlier, collisions with gaseous atoms randomize vectors of angular 
momenta to high degree. Due to efficiency of nuclear relaxation (LD99b) most 
of the grains do not rotate suprathermally in the absence of radiative torques. 
This explains why small grains for which radiative torques are not important 
are only marginally aligned. An alternative explanation of this fact is given by 
Mathis (1986) and is related to the preferential presence of ferromagnetic inclu- 
sions within large grains. However, it is claimed in Lazarian (2003) that data 
in Whittct et al. (2001) is consistent only with the change of the alignment ef- 
ficiency with the extinction rather than the grain size. More discussion of grain 
alignment can be found in the paper by Roberge (2004) in this volume. 

Structure of Grains 

While grain formation may favor loosely bound conglomerates of low fractal 
dimension, rapid rotation of grains destroys grain that are not sufficiently com- 
pact. The tensile stress calculated in Draine &: Lazarian (1998b) did not place 
severe constrains on small grain considered unless the grains consist of loosely 
connected parts. 

Implications of Translational Motion 

Shattering and CoagulaMon 

With the grain relative velocities known, we can make predictions for grain 
shattering and coagulation. For shattering, we adopt the Jones' et al. (1996) 
results, namely, the shattering threshold is 2.7km/s for silicate grain and 1.2km/s 
for carbonaceous grain. The critical sticking velocity were calculated in Chokshi 
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et al. (1993) (sec also Dominik & Tielens 1997).^ However, experimental work 
by Blum (2000) shows that the critical velocity is an order of magnitude larger 
than the theoretical calculation. Comparing these critical velocities with the 
velocity curve we obtained for various media, we can get the corresponding 
critical size for each of them (see Table 1). However, given the uncertainties 
with the parameters as discussed above, we are cautious about these numbers. 
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> 0.003 


> 0.001 


Coagulation size (/xm) 


< 0.01 < 0.02 


< 0.02 


< 0.05 


NA 


NA 



Table 1. The critical shattering and coagulation size in different 
medium. NA=not applicable. 



Correlation between turbulence and grain sizes 

Change in the intensity of turbulence should entail variations in grain sizes. 
The grain velocities are strongly dependent on the maximal velocity of tur- 
bulence V at the injection scale, which is highly uncertain. Thus the critical 
coagulation and shattering sizes would also depend on the amplitude of the 
turbulence accordingly. 

Composition of cosmic ray 

It has been shown that the composition of the galactic cosmic ray seems to 
be better correlated with volatility of elements (Ellison, Drury & Meyer 1997). 
The more refractory elements are systematically overabundant relative to the 
more volatile ones. This suggests that the material locked in grains must be 
accelerated more efficiently than gas-phase ions. This leads to the speculation 
of acceleration of grain erosion products in shocks (Epstein 1980; Cesarsky &: 
Bibring 1981; Bibring & Cesarsky 1981). If the sputtering happens upstreams, 
the sputtered products will be carried downstream, where they can be further 
accelerated to cosmic ray energies with higher efficiencies than gas-phase ions 
(Epstein 1980). The destruction by shocks may also be a prodigious source of 
PAHs, HACs and small grains (Jones, Tielens & Hollenbach 1996, Tielens et al. 
1994). The stochastic acceleration is more efficient for low speed grains and thus 
can serve as a preacceleration mechanism for shock acceleration. 

Vacuum cleaning and mechanical alignment 

Our results indicate that grains can get supersonic through interaction with 
fast modes. Grains moving supersonically can efficiently vacuum-clean heavy 
elements as suggested by observations (Wakkcr & Mathis 2000) . The supersonic 
grains can also be aligned (see a review by Lazarian 2003 and references therein). 
As pointed out earlier, the scattering is not efficient for slowly moving grains 
so that we may ignore the effect of scattering to the angular distribution of the 
grains. Since the acceleration of grains increases with the pitch angle of the 
grain, the supersonic grain motions will result in grain alignment with long axes 
perpendicular to the magnetic field. 



There are apparent misprints in the numerical coefhcient of Eq.(7) in Chokshi et al.(1993) and 
the power index of Young's modulus in Eq.(28) of Dominik & Tielens (1997). 
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Grain Segregation and Turbulent Mixing 

Our results are also relevant to grain segregation. Grains are the major 
carrier of heavy elements in the ISM. The issue of grain segregation may have 
significant influence on the ISM metallicity. Subjected to external forcing, e.g., 
due to radiation pressure^, grains gain size-dependent velocities with respect to 
gas. WDOl have considered the forces on dust grains exposed to anisotropic 
interstellar radiation fields. They included photoelectric emission, photodcsorp- 
tion as well as radiation pressure, and calculated the drift velocity for grains of 
different sizes. The velocities they got for silicate grains in the CNM range from 
O.lcm/s to lO^cm/s. Fig. 5a shows that the turbulence produces larger veloc- 
ity dispersions^'^. Those velocities are preferentially perpendicular to magnetic 
field, but in many cases the dispersion of velocities parallel^^ to magnetic field 
will be comparable or greater than the regular velocities above. 

More important is that if rcconncction is fast (see Lazarian & Vishniac 
1999), the mixing of grains over large scales is provided by turbulent diffusivity 
~ VL. Usually it was assumed that the magnetic fields strongly suppress the 
diffusion of charged species perpendicular to their directions. However, this 
assumption is questionable if we notice that motions perpendicular to the local 
magnetic field are hydrodynamic to high order as suggested by Cho, Lazarian 
& Vishniac (2002). In fact, recent work by Cho et al. (2003) shows that the 
diffusion processes in MHD turbulence are as efficient as in hydrodynamic case if 
the mean magnetic field is weak or moderately strong, i.e., Bq is ~ equipartition 
value. This means that from the theoretical perspective, we do expect that grains 
can be mixed by the MHD turbulence. Therefore we believe that the segregation 
of very small and large grains speculated in de Oliveira-Costa et al. (2002) is 
unlikely to happen for typical interstellar conditions. 

Interdependence of Rotational and Translational Motions 

Coupling of different types of motions is very important for understand- 
ing grain dynamics. As we discussed earlier, the coupling of vibrational and 
rotational motions results in flips that change dynamics of grains. 

Similarly translational and rotational motions are interdependent. For in- 
stance, frequent flips average out not only uncompensated torques, but also the 
uncompensated force. As we mentioned earlier, radiative torques can rescue 
grains from thermal trapping and therefore we show in Fig 5 that grains larger 
than ~ 5 X 10~^ cm experience uncompensated forces due to H2 formation. 

Fast rotation makes atoms less susceptible to mechanical alignment through 
gaseous collisions. However, as grains undergo crossovers, they can still be 



Even in the absence of radiation pressure grains can move along magnetic field lines due to 
the uncompensated forces, e.g. due to an unequal rmmbcr of active sites of H2 formation (see 
P79). Those forces would be mitigated in molecular clouds, which would induce inflow of dust 
into molecular cloud. The latter would affect metallicity of the newborn stars. 

^''Our calculation show that for the chosen set of parameters the effects of H2 thrust are also 
limited. 

^^This dispersion stems from both the fact that the transpositions of matter by fast modes 
are not exactly perpendicular to magnetic field (see plot in Lazarian & Yan 2002b) and due 
to randomization of directions of grain velocities by magnetized turbulence (Yan & Lazarian 
2003). 
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aligned via processes that were termed in Lazarian (1995a) crossover and cross 
sectional alignment (see also Lazarian & Efroimsky 1996; Efroimsky 2002). 

Streaming grains get rotation rates that are faster than that of their ther- 
mally rotating counterparts. This may allow thermally trapped grain get "un- 
trapped" and therefore acquire suprathermal rotation through uncompensated 
Purcell's torques (P79). If grains are "helical" in terms of collisions with gaseous 
atoms the effect of grain spinning up and getting untrapped should be even move 
pronounced. 

In addition, some characteristics of dust grains are affected both by rota- 
tional and translational motions. For instance, both fast rotation and collisions 
limit the structure of grains and prevent very loose aggregates to dominate the 
extinction. 

7. Summary 

1. Interstellar grains move relatively to each other and this results in collisions. 
The outcome of those may be shattering or coagulation depending on the relative 
velocities of grains. As grains move relatively to gas they adsorb atoms from the 
gas and this affects dust chemistry. In addition grains get aligned as a result of 
gas- dust streaming. Grain rotation affects both grain alignment and grain size 
distribution. Fast rotating grains are not subjected to randomization via gas 
collisions. Fast rotation destroys loosely connected grains. 

2. Grains do not always rotate about their axis of major inertia. Ther- 
mal fluctuations cause vibrations of grain angular momentum about the axis of 
major inertia and occasional flips. The rate of flips depends on the coupling of 
rotational and vibrational degrees of freedom that is determined by the rate of 
dissipative processes in the wobbling grain. Frequent flips result in grains being 
thermally trapped, i.e. rotating with thermal velocities in spite of the presence 
of suprathermal torques. 

3. The rate of grain rotation is determined by collisions with neutrals, 
ions, photons, plasma effects, emission of radiation, nascent II2 molecules etc. 
It is also influenced by the grain flipping. Depending on the interplay of these 
processes, grains can rotate both at a rate that is larger and smaller than the 
rate of the thermal Brownian motion. 

4. AlthoTigh rotation is classical even for the smallest grains, a number 
of STibtlc quantum mechanical effects determine grain dynamics. For instance, 
plasma drag arising from the interaction of the grain dipole moment with imping- 
ing ions is limited by the quantum nature of the interactions. Nuclear internal 
relaxation is a quantum effect arising from reorienting of spins of nuclei within a 
wobbling grain. Paradoxically, in terms of coupling of rotational and vibrational 
degrees of freedom nuclear spins can be much more efficient than their electron 
counterparts. Resonance paramagnetic relaxation is yet another quantum effect 
that alleviates alignment of tiny grains. 

5. Grains can be accelerated by various mechanisms. Formation of H2, 
variations of the accommodation coefficient, photoelectric emission and pho- 
todesorption, radiation pressure are the processes that can induce grain motion. 
However, ubiquitous interstellar turbulence is probably the major driver of rela- 
tive grain-gas and grain-grain motions. As the turbulence is usually magnetized 



18 



Lazarian &; Yan 



and grains are usually charged the magnetic nature of turbulence is essential. 
Magnetic fields, on one hand, limit the relative gas-grain motions in the direc- 
tion perpendicular to the local magnetic field. On the other hand, interaction 
of charged grains with magnetic turbulence results in resonance acceleration of 
dust. Supersonic grain velocities are attainable as a result of such an interaction. 
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